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ABSTRACT 

A supermassive black hole in the nucleus of a galaxy will be revealed when a 
star passes close enough to be torn apart by tidal forces and a flare of radiation 
is emitted by the stream of stellar debris that plunges into the black hole. Since 
common active galactic nuclei have accreting black holes that can also produce 
flares, a convincing demonstration that a stellar tidal disruption has occurred 
generally begins with a "normal" galaxy that has no evidence of prior nuclear 
activity. Here we report a luminous UV flare from an elliptical galaxy at z — 0.37 
in the Groth field of the GALEX Deep Imaging Survey that has no evidence of 
a Seyfert nucleus from optical spectroscopy and X-ray imaging obtained during 
the flare. Multiwavelength data collected at the time of the event, and for 2 years 
following, allow us to constrain, for the first time, the spectral energy distribution 
of a candidate tidal disruption flare from optical through X-rays. The luminosity 
and temperature of the radiation and the decay curve of the flare are in excellent 
agreement with theoretical predictions for the tidal disruption of a star, and 
provide the strongest empirical evidence for a stellar disruption event to date. 
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Introduction 



Dynamical studies of nearby galaxi es show that all galax ies with a bulge host a su- 
permassive black hole in their nucleus (jMagorrian et al. 19981) . Stellar dynamical models 
(IMagorrian &: Tremaine 19991 ; IWang fc Merritt 2004 ) predict that once every 10 4 -10 5 yr the 
orbit of a star in the nucleus of a galaxy will pass within the tidal disruption radius of its cen- 
tral supermassive black hole, Rt ~ ^(Mbh/M*) 1 / 3 , an d the star will be disrupted. A stellar 
disruption results in a flare of electromagnetic radiation (IFrank fc Rees 19761 ; iLidskii fc Ozernoi 19791 ) 
as the bound fraction of the stellar debris, ^ 0.5M* (Rees 1988; Ayal et al. 2000), falls back 
onto the black hole and is accreted. There is a critical black hole mass above which Rt is 
smaller than the S chwarzschild radius (R s ), and the star is swallowed whole without dis- 
ruption (IHills 19751 ). For a solar-type star, M cr u ~ 1O 8 M . The luminosity and decay of a 
tidal disruption flare is dependent on the mass and spin of the central black hole, and can 
be used to directly probe dormant black holes in galaxies for which the sphere of influence of 
the black hole (R sp h = GMbh/c*) is unresolved, and a dynamical measurement of the black 
hole mass is not possible. 

The start of the flare, to, occurs when the most tightly bound debris first returns to 
pericenter (R p ) after the star is disrupted at time to, where 

(t -t D ) = 0Alk- 3 /\R p /R T )\RjR & f/\MjM & )-\M Bll /10 G M Q )^ yr 

(Li et al. 2002), and k — 1 for a star with no spin, and k = 3 for a star that is spun 
up by tidal interactions to near break-up on disruption. After to, the return rate of the 
bound debris to pericente r declines as function of time as M{t) oc [(t — t D )/(t — t D )]~ 5 / 3 
(jEvans &: Kochanek 19891 ). which determines the characteristic power-law decay of the lu- 
minosity of the flare over the following months and years. When the mass accretion rate is 
super-Eddington, 

M Edd = L Edd /(ec 2 ) = O.O25(M BH /lO 6 M )(e/O.l)- 1 M yr- 1 , 

where e is the efficiency of converting mass into radiation, the debris is accreted in a 
thick disk that radiates close to the Eddington luminosity (jUlmer 19991 ). LEdd = 1.3 x 
10 44 (Mbh/10 6 M q ) ergs s -1 . The spectrum of a tidal disruption event can thus be charac- 
terized by the blackbody temperature of a thick disk radiating at L^m at radii ranging from 
Rt to i? p , T e fj « (L F/M /4 :iiaR 2 ) 1 / 4 ', w hich for a solar-type star is ~ (2 — 5) x 10 5 K, and 
peaks in the extreme-UV (jUlmer 19991 ) . 



The most unambiguous cases for a stellar disruption occur from host galaxies with no 
evidence of an active galactic nucleus (AGN) for which an upward fluctuation in the ac- 
cretion rate could also explain a luminous UV/X-ray flare. A UV flare from the nucleus 
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of the elliptical galax y NGC 4552 was pr oposed to be the result of the tidal stripping of 
a stellar atmosphere (jRenzini et al. 19951 ); however the possible presence of a persistent, 
low-luminosity AGN detected in hard X-rays (Xu et al. 2005) makes this interpretation 
uncertain. The ROSAT All-Sky survey conducted in 1990-1991 sampled hundreds of thou- 
sands of galaxies in the soft X-ray band, and detected luminous (10 42 — 10 44 ergs s -1 ), soft 
[Tbb = (6 — 12) x 10 5 K] X-ray flares from several galaxies with no previous eviden ce for 
AGN activity, and with a flare rate of ~ 1 x 10~ 5 yr" 1 per galaxy (IDonley et al. 20021 ). that 
is consistent with the theoretical stellar disruption rate. A decade later, follow-up Chan- 
dra and XMM-Newton observations of three of the galaxies demonstrated that they had 
faded by factors of 240 — 6000, consistent with the (t — t D )~ 5 ^ 3 decay of a tidal disruption 
flare (Komossa et al. 2004; Halpern et al. 2004). Follow-up Hubble Space Telescope (HST) 
Space Telescope Imaging Spectrograph narrow-slit spectroscopy confirmed two of the galax- 
i es as inactive, qu alifying them as the most convincing hosts of a tidal disruption event 
( jGezari et al. 20031 ) . The ROSAT flare with the best sampled light curve was successfully 
modeled as the tidal disruption of a brown dwarf or planet (Li et al. 200 2), although its hos t 
galaxy was subsequently found to have a low-luminosity Seyfert nucleus (jGezari et al. 20031 ) . 



2. Observations 

We initiated a program to take advantage of the UV sensitivity, large volume, and 
temporal sampling of the Galaxy Evolution Explorer (GALEX) Deep Imaging Survey (DIS) 
to search for stellar disruptions in the nuclei of galaxies over a large range of redshifts. The 
DIS covers 80 deg 2 of sky in the far-ultraviolet (FUV; A = 1344 — 1786A) and near-ultraviolet 
(NUV; A = 1771 — 283lA) with a total exposure time of 30-150 ks, that is accumulated in 
~ 1.5 ks eclipses (when the satellite's 98.6 minute orbit is in the shadow of the Earth). 
Due to target visibility and mission planning constraints, some DIS fields are observed over 
a baseline of 2-4 years to complete the total exposure time. This large range in cadence 
of the observations allows us to probe variability on timescales from hours to years. Here 
we present a UV flare discovered in the GALEX 1.2 deg 2 Groth field at a position of right 
ascension 14 h 19 m 29. s 8, declination +52°52'6". We use multiwavelengt h survey data from the 



All- Wavelength Extended Groth Strip International Survey (AEGIS) (IDavis et al. 20061 ) and 
the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) to identify the host of the 
flare, and study the broadband properties of the flare and its subsequent decay. Throughout 
this Letter, calculations are made using Wilkinson Microwave Anisotropy Probe cosmological 



parameters (IBennett et al. 20031 ): Hq = 71 km s 1 Mpc 1 , Qm — 0.27, and Q\ = 0.73. 



The UV flare reported here was discovered as a result of its large amplitude of variability 



-4- 



during 4 years of DIS observations. Figure [T] shows the discovery images in the NUV and 
FUV of the source, with the eclipses co-added into yearly epochs, and Figure [2] shows its 
detailed NUV+FUV light curve. The source is undetected in the NUV and FUV in 2003 
June 21 - 29, and then appears in 2004 March 25 - June 24, indicating an amplitude of 
variability of Am ^ 2 mag. The source then decays monotonically by ~ 2.0 mag over 2 
years to almost below the detection threshold in 2006 March 5-7. Archival AEGIS 50 ks 
Chandra (0.3 - 10 keV) observations on 2005 April 6 and 7 detected an extremely soft X-ray 
source that appeared on the second day of the observations, with 10 photons with energies 
between 0.3 and 0.8 keV, indicating a spectral slope, Je oc E~ t , of T = 7 ± 2 (reduced 
X 2 = 0.92), for a column density of neutral hydrogen fixed to the Galactic value. No source 
was detected in the following 50 ks Chandra observations on 2005 September 20 and 23. 



3. Interpretation 

The UV flare and variable soft X-ray source are coincident with an elliptical galaxy 
at z = 0.3698 (luminosity distance, d L = 1970 Mpc). Figure [3] shows the AEGIS HST 
ACS optical image with the UV and X-ray positions overplotted, and the AEGIS Keck 
DEIMOS optical spectrum from the DEEP2 survey of the host galaxy taken on 2005 March 
9. The spectrum has strong stellar absorption features typical of an early-type galaxy. We 
marginally detect faint narrow [O III] A5007 line emission with L([0 III]) = 8 ± 4 x 10 39 
ergs s -1 . The non-detection of a hard X-ray source by Chandra places an upper limit to 
the 3-20 keV luminosity of a power-law AGN spectrum in the nucleus of Lx ^ 1 x 10 42 
ergs s~ l . The correla tion between L([0 III]) and Lx(3-20 keV) observed for local AGNs 



( jHeckman et al. 20051 ) predicts L([0 III]) < 10 39 ergs s^ 1 , which suggests that the possible 
[O III] line emission has too high an L([0 III])/Lx ratio to be associated with AGN activity, 
and may be from low levels of star formation. We obtained a follow-up spectrum after 
the flare discovery on 2006 June 1 with the Double Spectrograph on the Palomar 200 inch 
telescope, and again detected no strong emission lines. CFHTLS variability monitoring data 
from 2005 January to June show no detection of a variable optical source in the galaxy, with 
up to 5 observations per month with a sensitivity of mu m = 25. 

The lack of variable optical emission during the UV flare puts an important constraint 
on its nature. Supernovae (SNe) have light curves that extend over months, similar to 
the flare observed here, but their emission peaks at optical wavelen gths, and is intrinsi- 



cally faint shortward of ~ 3000 A ( Panagia 20031 ; iBrown et al. 2005 ). Gamma ray bursts 



(GRBs) have power-law decays, but again, their emission is brighter at optical wavelengths 



( IMeszaros fc Rees 19971 ). The characteristic power-law continuum of AGN emission ex- 
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tends from optical to extreme- UV as /„ oc v a , where 0.43 < a < 2 (IZheng et al. 1997 



Dietrich et al. 20021 ). If the UV flare were an upward fluctuation of AGN continuum emis- 
sion, a power-law spectrum with a > 0.43 would produce an optical source with m < 23 
mag during the C FHTLS observations and power a broad H/3 line with a luminosity of ^ 



1 x 10 41 ergs s -1 ([Dietrich et al. 20021 ) at the time of the Keck spectrum, both of which are 
definitively ruled out by the observations. In summary, the flare is most probably not due 
to AGN activity because it displays (1) no optical variability, (2) no broad emission lines, 
(3) a soft X-ray spectrum with no hard X-ray component, and (4) an early- type galaxy 
spectrum with possible narrow [O III] emission with a high [O III]/Lx ratio characteristic of 
star formation. In addition, a SN or GRB origin is unlikely because of the high UV/optical 
ratio. 

The spectral energy distribution (SED) of the flare from optical to X-rays on 2005 April 
7 is shown in Figure HI A blackbody fit to the SED is strongly constrained by the soft X-ray 
flux density, which is best fit for a rest-frame = 4.9 x 10 5 K, shown with a solid line 
in Figure HI which for L Ve= (i +z ) Vo = f Vo (ATtd 2 L )/{l + z), corresponds to Rbb = 6.2 x 10 12 cm 
and L bo i = 1.2 x 10 45 ergs s _1 . The steep Wein's tail of the blackbody curve is consistent 
with the extremely soft spectral slope of the X-ray source. However, since the X-ray flux 
varies on a short timescale, and only appears above the detection threshold for part of the 
observations, we consider this fit an upper limit to for the duration of the flare. The 
lowest Tbt, that is consistent with the UV flux densities and the optical upper limits is ~ 1 
x 10 5 K, shown with a dashed line in Figure HI which corresponds to Rbb = 1.6 x 10 13 cm 
and Lboi = 1.3 x 10 43 ergs s _1 . This range of temperatures and luminosities is in excellent 
agreement with the theoretical predictions for a tidal disruption flare. 



4. Discussion 

The characteristic tidal disruption (t — tf))~ 5 ^ 3 power-law decay fit to the UV light curve 
yields a best-fit time of disruption of t£)=2003.3 ± 0.2. If we allow the power-law index of 
the decay, (t — to) ~ n , to vary as well as to, we get a comparable fit, but with much less 
constrained values for to = 2002 ± 2 and n = 3 ± 2. The GALEX observations do not cover 
the rise of the UV flux to its peak value, however we can constrain the time of the peak 
(to) to be between the time of the GALEX detection beginning on 2004 March 25 and the 
latest date of the non-detection on 2003 June 29, (t„ - t D ) = 0.2(±0.2) - 0.9(±0.2) yr. We 
estimate the total energy released during the flare by integrating the best-fit (t — £d)~ 5 / 3 
light curve from the time of the peak of the flare to infinity, E to t = L(t)dt. If on 2005 
April 6-7, 1.3 x 10 43 ergs s" 1 < L M < 1.2 x 10 45 ergs s" 1 , and (t„ -°t D ) < 1.1 yr, this 
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yields a minimum total energy released of 1.8 x 10 51 ergs £ E tot £ 1.7 x 10 53 ergs, and a 
minimum mass accreted, M acc = E tot /(ec 2 ), of O.OlM (e/O.l) _1 < M acc < O.9Af (e/O.l) _1 , 
assuming T bb = (1 — 5) x 10 5 K during the flare. 

For a solar-type star, the limits on (t — t D ) imply k 3 3 x 1O 6 M £ M BH £ k 3 7 x 1O 7 M , 
where k depends on the spin-up of the star, with a hard upper limit on the black hole 
mass of M BH < M crit = 1.15 x lO 8 (r 3 /m) 1/2 M , where r = R*/R & , and m = M*/M Q , 
in order to satisfy the condition that Rp > R s , and a stellar disruption can occur. The 
radius of the flare emission also places constraints on Mbh- The debris disk should form 
within the tidal disruption radius of the central black hole. The blackbody fits to the flare 
SED have R bb ^ 16 x 10 13 cm which, if R bb ~ Rt, places a lower limit on the black hole 
mass of Mbh ~ 12 x 1O 7 M . The radius of the emission should also be greater than the 
minimum stable particle orbit for the black hole (R ms ), which ranges from R ms = 6R g for 
a black hole with no spin, down to R ms = R g for a maximally spinning black hole, where 
R g = GM B n/c 2 . For R bb > R ms , this requires that M B h < 11 x 1O 8 M , and is consistent 
with Mbh < M crit for a solar-type star. We can also use the strong co rrelation between blac k 



hole mass and stellar velocity dispersion established for local galaxies (iTremaine et al. 20021 ). 
log(M BH /M ) ~ 8.13 + 4.021og(a e /200 km s" 1 ), to estimate M BH - The stellar velocity 
dispersion within the half-light radius (R e ) of the bulge of the galaxy, a e = 120 ± 10 km 
s , yields a mass of 2±l x lO 7 Af , where the 1 a error includes the intrisic scatter of the 
correlation of ~ 0.3 dex. 

To summarize, the black hole mass estimated from the Mbh — 0" relation (Mbh = 
2~tl x 1O 7 M ) is consistent with the constraints on the black hole mass from the blackbody 
spectrum such that R ms < R bb <, R T (1.2 x 1O 7 M < M B h < 11 x 1O 8 M ), with the limits 
derived from the time delay between the time of the disruption and the peak of the flare 
(k 3 3 x 1O 6 M < M BH ~ k 3 7x 1O 7 M ), and with the condition that M BH < M crit = 1.15 x 10 8 
for a solar-type star to be disrupted. It is encouraging that these various independent black 
hole mass estimates are in agreement with each other as long as k < 3, and the star was not 
spun-up to near break-up on disruption. The GALEX DIS has proven to be a promising 
survey for detecting stellar disruption events so far. With a larger sample of tidal disruption 
flares, we can begin to measure the distribution of masses and spins of black holes in the 
nuclei of normal galaxies as a function of redshift, unbiased by the minority of galaxies that 
are AGNs. 
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Fig. 1. — GALEX images of the flaring source over 4 years of DIS eclipses coadded into 
yearly epochs of NUV (t exp = 16, 1.5, 110, and 29 ks) and FUV (t exp = 16, 1.5, 41, and 29 ks) 
observations, labeled with the mean date of each coadd. Blue circles show the 9" diameter 
aperture in the NUV, and the 12" diameter aperture in the FUV used for photometry. No 
source is detected in NUV and FUV in the 2003.48 epoch. A smaller aperture was used 
in the NUV to avoid contamination by a pair of foreground interacting galaxies detected in 
the NUV just South of the source. The counts from the background in the aperture were 
determined from the pipeline background image, and an aperture-correction factor of 80% 
and 85% was applied to the NUV and FUV aperture fluxes respectively, to recover the total 
energy enclosed in a point source. 
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Fig. 2.— GALEX light curve in the FUV (A e// = 1528A) and NUV (A e// = 2271 A). Error 
bars show 68% confidence intervals and arrows show 95% confidence upper limits based on 
Bayesian statistics. For the later observations in 2005 the FUV detector was temporarily 
not operational. The least-squares fit to the NUV light curve with a (t — £d) ~ 5 ^ 3 decay is 
shown with a thick solid line, yielding a best fit time of disruption to = 2003.3, with a 1 
a error of 0.2 yr. The least-squares fit to the NUV light curve with the power-law index, 
(t — t D )~ n , and t D allowed to vary is shown with a dashed line, with t D = 2002 ± 2 and 
n = 3 ± 2. Dotted lines with an X indicate times of two sets of 100 ks Chandra 0.3 - 10 keV 
observations. An extremely soft X-ray source was detected during the first set observations, 
and no source was detected in the second set of observations. CFHTLS monitoring data 
from 2005 January to June in g (A e// = 4763A),r (A e// = 6174A),i (A e// = 7619A), and 
z (A e // = 8847A) with mu m = 25, detect no variable optical source during the UV flare. 
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Fig. 3. — Image: AEGIS HST ACS /-band image (A e // = 8140A) of the galaxy host with 
the GALEX 1" position error circle {blue) of the UV flare and the Chandra 0."7 position 
error circle (red) of the soft X-ray source, with a 1" systematic astrometry correction applied 
to the Chandra position measured from the shift of the Chandra astrometry with respect to 
ACS for a nearby star. The morphology of the galax y is determine d from a two-dimensional 
bulge/disk decomposition of the ACS I-band image (ISimard 19981 ). which results in a bulge 
fraction of 0.72, and a half-light radius of the bulge, R e = 0."43. Spectrum: AEGIS Keck 
DEIMOS spectrum from the DEEP2 survey of the galaxy host obtained on 2005 March 
9, smoothed by 6 pixels (~ 2 A), and not corrected for the instrumental response. The 
spectrum shows strong stellar absorption lines typical of an early-type galaxy, with a marginal 
detection of narr ow [O HI] A5007 line emi ssion. The best-fitting Bruzual- Chariot early-type 
galaxy template (IBruzual fc Chariot 20031 ) scaled to the continuum of the spectrum is shown 
in green. The luminosity-wei ghted line of sight velocity dispersion measured using stellar 
templates JlVforan et al. 2005^ within the 1" slit width of DEIMOS is 120 ± 10 km/s, which 
corresponds to an ape rture of 1.2R e , for which the correction to an aperture of R e is < 5% 
jGebhardt et al. 200oh . 
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Fig. 4. — Spectral energy distribution from optical to X-rays of the flare on 2005 April 7, 
with the best blackbody fit with T&& = 4.9 xlO 5 K shown with a solid line. The spectral slope 
of the extremely soft X-ray detection, /„ oc E~ r , of T = 7 ± 2 is also shown, corresponding 
to an flux density at 0.5 keV of (1.3 ±0.8) x 10~ 31 ergs s _1 cm -2 Hz -1 , a factor of ~ 4 
above the detection limit of Chandra (dotted line). Since the X-ray flux appears above the 
detection threshold of Chandra for only part of the observations, the blackbody fit gives an 
upper limit to during the flare. Dashed line shows the lowest temperature blackbody 
consistent with the UV flux densities and optical upper limits, with T bb = lxlO 5 K. 



